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Abstract 

We report formulas for the joint moments of the determinantal products (det p)'^(det p^'^Y — 
0, 1, 2, . . . , A^; K = 1, . . . , 12) of Hilbert-Schmidt (HS) probability distributions over the generic two- 
rebit and two-qubit density matrices p {k = 1, ... ,4). Here PT denotes the partial transposition 
operation of quantum-information-theoretic central importance. Each formula is the product of the 
expression for the HS moments of (det p)^, k = 0,1,2, ... , A^~special cases of results of Cappellini, 
Sommers and Zyczkowski (Phys. Rev. A 74, 062322 (1996))-and an adjustment factor. The 
factor is a biproper rational function, with its numerators and denominators both being 3K-degree 
polynomials in k. We infer the structure that the denominators follow for arbitrary k in both 
the two-rebit and two-qubit cases, and the six leading-order coefficients of k of the numerators in 
the two-rebit scenario. We also commence an analogous investigation of generic rebit-retrit and 
qubit-qutrit systems. This research was motivated, in part, by the objective of using the computed 
moments to well reconstruct the HS probabilities over the determinant of p and of its partial 
transpose, and to ascertain-at least to high accuracy-the associated (separability) probabilities of 
"philosophical, practical and physical" interest that (detp^"^) > 0. 
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We begin our investigation into certain statistical aspects of the "geometry of quantum 
states" PQ 12] by noting the two following special cases-which we will extend below-of the 
general formulas [3][eq. (3.2)]: 



and 

A; = 0, 1, 2, . . . The bracket notation () is employed by us to denote expected value, while p 
indicates a generic (symmetric) two-rebit or generic (Hermitian) two-qubit (4 x 4) density 
matrix. The expectation is taken with respect to the probability distribution determined 
by the Hilbert-Schmidt/Euclidean/fiat metric on either the 9-dimensional space of generic 
two-rebit or 15-dimensional space of generic two-qubit systems [Hll]. 

We report below sixteen (twelve two-rebit and four two-qubit) non-trivial extensions of 
these formulas, involving now in addition to |p|, the quantum-theoretically important deter- 
minant Ip'^^I of the partial transpose of p. (The nonnegativity of Ip'^^j-by the celebrated 
Peres-Horodeccy results [SH3|-constitutes a necessary and sufficient condition for separabil- 
ity/disentanglement, when p is either a4x4or6x6 density matrix.) At this point of our 
presentation, we note that three of these extensions are expressible-incorporating on their 
right-hand sides the two formulas above-as 



,fc|PT. (fc-l)(fc(2fc + ll) + 16) 

_ k{k{k{k{4:k{k + 12) + 203) + 368) + 709) + 2! 
2-rebit/HS - io24(A; + 3)(A; + 4)(4A; + U){Ak + 13){Ak + 15){4k + 17) ^1^1 ^2-refeit/H5 



32(A; + 3)(4A; + 11)(4A; + 13) 
, pr,2, fc(fc(fc(fc(4A;(A; + 12) + 203) + 368) + 709) + 2940) +4860 ,, .f^. 

\P\ \P I )2-rebit/HS = mo^^ i o\n i A\rA1 i i T o\ r A 1 i T^\rA1 I 1^7^ \IPl / 



(4) 



and 



/\^\k\PT\\ _ k{k{k + 6) 1) 42 fr\ 

{\P\\P \)2-,ubu/HS - g^2A; + 9)(4fc + i7)^4k + 19) ^^-.^..^/h^- (5) 

These three formulas were, first, established by "brute force" computation-that is calcu- 
lating the first {k = 0, 1, 2, . . . , 15 or so) instances, then employing the Mathematica com- 
mand FindSequenceFunction, and verifying any formulas generated on still higher values of 
k. {Initially, although we had the specific values of {\p\'^\p^'^\^)2~rebit/HS for = 0, . . . , 13, 
and similarly for {\p\''\p^'^\^)2-rebit/HSy "we were not able to determine, in the same manner, 
encompassing expressions for them.) 



2 



As a special case (fc = 1) of ([3]), we obtain the rather remarkable moment result, zero, 
already reported in |8j. The immediate interpretation of this finding is that for the generic 
two-rebit systems, the two determinants \p\ and \p^'^\ comprise a pair of nine-dimensional 
orthogonal polynomials [914TT] with respect to Hilbert-Schmidt measure. (C. Dunkl has 
kindly pointed out that orthogonality here does not imply zero correlation.) In addition 
to this first {k = 1) HS zero-moment of the ("equally-mixed") product variable |p||p^"^| in 
the two-rebit case, we had been able to compute its higher-order moments, k = 2, ... ,6. 
(The results for k = 2 can be obtained by direct application of (|4]). The feasible range of 
the variable is G [~ 110592 ' 2^]^^^^ lower bound of which— j^;^^ = — 2^^^3^^ we 

determined by analyzing a general convex combination of a Bell state and the fully-mixed 
state.) 

These five further moments of IpUp"^"^!, k = 2, ... ,6, are all rational numbers. T If we 
take the ratios of these first six moments of |p||p'^"'"| to the first six even moments given by 
([1]), that is the values {\p\'^'^)2-rebit/HS, k = 1, ... ,6, we obtain the rather succinct sequence, 

{\pn2-reMt/HS ^ ' 54' 55' 175' 833' 3724^ ^ ^ 

^ {0, 1.425926, 0.4363636, 1.194286, 0.7178872, 1.055048}. 

(As to the two-qubit counterpart of this sequence, we had so far only been able to compute 
its very first term-turning out, quite remarkably, to be the negative value — |.) 

Since these ratios ^ are so comparatively simple, it suggested to us that we might be 
more able to progress in a series of analyses [T2HT9] (mainly devoted to the determination of 
separability probabilities), by making our initial goal the computation of these ratios for still 
higher-order moments-rat her than the direct computation of the very small values, having 
lengthy multi-digit denominators, of the moments {{\p\\p^'^\)'')2-rebit/HS themselves. (In 
[8][eqs, (33)-(41)], we were able to report and analyze the first nine moments (|p^^|'^)-the 
first two of which, — ^ and „ . . . , can be obtained by directly setting = in (p3|) and (0), 



858 2489344' _ _ 

respectively. However, to this point, we have not found any associated similarly compact 
sequences of moment ratios, as above.) 

we display the sequence of ratios — (|pP'°)2 ^t't/gs^^ ' 
k = 1, . . . , 100, the first six members of which have been exactly calculated, 
as noted above, and the rest through extended-precision (60-digit) numerical com- 
putations. (Simply as an indicator of accuracy of these computations, the nu- 
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FIG. 1: The six- member two-rebit HS exact moment ratio sequence ([6]), supplemented by its 
numerical continuation, using extended-precision (60-digit) arithmetic. Three hundred and twenty 
million random density matrices were employed. 

merical estimates yielded by this procedure of the six-member exact sequence ([G]) 
were {-0.0002052822,1.426286,0.4359643,1.194784,0.7175908,1.055326}.) Similarly, in 
Fig. |2| we display the (quite differently-behaving) two-qubit sequence of moment ratios 
^^'ll^Ifcr^'"""'"^"'"' , = 1, . . . , 100, only the first member of which, -|, we were initially able 

\\P\ /2 — qubit/HS ^ 

to exactly compute, and the remaining ninety-nine, numerically, using extended-precision. 

If we are, at some point, in the course of these extended analyses, able to develop formulas 
explaining the jull sequences of ratios in the two-rebit and two-qubit cases, we should be 
able to reconstruct the Hilbert- Schmidt univariate probability distributions over the product 
variable |p||p^"^| (cf. [3j[Figs. 2-4]). From such reconstructed distributions, HS separability 
probabilities should be determinable to high accuracy. 

For the further edification of the reader, we present in Fig. |3] a contour plot of the joint 
Hilbert-Schmidt (bivariate) probability distribution of \p\ and \p^'^\ in the two-rebit case, 
and in Fig. |4| its two-qubit analogue. (A colorized grayscale output is employed, in which 
larger values appear lighter.) In Fig. [sjis displayed the difference obtained by subtracting the 
second (two-qubit) distribution from the first (two-rebit) distribution. (The black curves in 
all three contour plots appear to be attempts by Mathematica to establish the nonzero- zero 
probability boundaries-which, it would, of course, be of interest to explicitly determine, if 
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FIG. 2: The two-qubit analogue of the two-rebit sequence depicted in Fig. [T] with only the first 



member 



having initially been exactly known, and the next ninety-nine computed numerically, 



using extended-precision (60-digit) arithmetic. Twenty-four million random density matrices were 
employed. 
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FIG. 3: Contour plot of the joint Hilbert-Schmidt probability distribution of |p| (horizontal axis) 
and \p^'^\ in the two-rebit case. Larger values appear lighter. The variable ranges are |p| G [0, 



and \p^'^\ G 



16' 256 J 



One billion random density matrices were employed. 
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FIG. 4: Contour plot of the joint Hilbert-Schmidt probability distribution of \p\ (horizontal axis) 
and \p^'^\ in the two-qubit case. Six hundred million random density matrices were employed. 

possible-of the joint domain oi\p\ and \p^'^\.) 

These last three figures are based on Hibert-Schmidt sampling (utilizing Ginibre ensem- 
bles [3]) of random density matrices, using 10,000 = 100^ bins. In regard to the two-qubit 
plot, K. Zyzckowski informally wrote: "A high peak in the upper corner means that: a) a 
majority of the entangled states is 'little entangled' (small det{p^)) or rather, they are 'close' 
to the boundary of the set, so one eigenvalue is close to zero, and the determinant is small; 
b) as det{p) is also small, it means that these entangled states live close to the boundary of 
the set of all states (at least one eigenvalue is very small), but this is very much consistent 
with the observation that the center of the convex body of the 2-qubit states is separable 
(so entangled states have to live 'close' to the boundary). Similar reasoning has to hold in 
the real case as well." 

At a later point in our investigation, we realized that we might make further progress- 
despite limitations on the number of moments we could explicitly compute-by exploiting 
the evident pattern followed by our newly-found formulas ^ and (|4])-in particular, the 
structure in their denominators. This encouragingly proved to be the case, as we were able 
to establish that 
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FIG. 5: Difference obtained by subtracting the two-qubit HS probability distribution in Fig. [4]from 
the two-rebit probabihty distribution in Fig. [3j Darker colors indicate more negative values. 

where 

A3 = 8A;^+180A;®+1674A;^+8559A;^+29493A;^+84291A;^+136801A;^-401334A;2-2516616A;-3612816 



(8) 



and 



B3 = 32768(A; + 3){k + 4) {k + 5) (4A; + 11) (4A; + 13) (4A; + 15) (4A; + 17) (4A; + 19) (4A; + 21) . (9) 



So, it is now rather evident that we can write for general non-negative integer k, 

A 

\\P\ IP I )2-rebit/HS — -f^\\P\ l2-rebit/HSi 



(10) 



where both the numerator and the denominator are 3K-degree polynomials (thus, 
forming a "biproper rational function" [20]) in k (the leading coefficient of being 2'^), 
and 



= 128'^(fc + 3)« 2k 
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(11) 



2k 



where the Pochhammer symbol (x)„ = ^^^^^^ = x{x + l) . . .(x + n — 1) is employed. Further 
still, moving upward to the next level (k = 4), we have determined that 



\P\ IP I 12-rebit/HS — 75~\IPI 1 2-rebit/ H S i 



(12) 



where 

= 16k^^ + 576A;" + 9112^^° + 84496A;^ + 525681 A;^ + 2389416A;^ + 7805462A;^ + 13904508A;^+ 

(13) 

+6212189A;^ + 166748972A;^ + 1636873812/^^ + 5496485760A; + 6610161600, 



and -B4 is given by (11) with k = 4. The real part of one of the roots of A4 is 2.999905, 
suggesting to us some possible interesting asymptotic behavior of the roots of these numer- 
ators, K, —7- 00. In our previous related study [8] [sec. II. B. 2], we were also able to discern the 
general structure that the denominators of certain "intermediate [rational] functions" used 
in computing the (univariate) moments of {p^'^\'^)2-rebit/HS, k = 1, . . . ,9 followed. 

From our four new two-rebit moment results ([3]), (|4]), ([T]) and (12), we see that the 
constant terms in the 3/t-degree numerator are —16,4860,-3612816 and 6610161600 
for K = 1,2,3,4. Since we had previously computed [8][eqs, (33)-(41)] the moments of 
{\p^'^\'^)2-rebit/HSy K = 1, . . . , 9, wc are also able to determine the next five members of this 
sequence { — 16,4860,-3612816,6610161600}. However, no general rule for this sequence, 
which would directly allow us to obtain a formula for {\p^'^\'^)2-rebit/HS, has, to our dis- 
appointment, yet emerged for them. (With such a rule, we could address the separability 
probability question through the reconstruction of a univariate probability distribution.) 

A simple algebraic exercise involving ([T]) shows that if we multiply the conversion factor 

^ ^ v/^r(2/. + 4)r(8/. + io) 

8r (2k + I) r(4ft: + 2)r(4/t + 10) ^ ^ 

by the rational function factors ^ found above that applied to {\p\'^)2-rebit/HS yield 
{{\p\\p^'^\)'^)2-rebit/HS, we obtain the ^-member of the sequence of moment ratios ([6|. Since 
the members of this sequence appear (Fig. [T]) to asymptotically approach 1 (our numerical 
estimate for the 100-th term is 1.001542), it would seem that the conversion factor c and ^ 
asymptotically approach one another. 

Certainly, it would be of interest to conduct analyses parallel to those reported above 
for metrics of quantum-information-theoretic interest other than the Hilbert- Schmidt, such 
as the Bures (minimal monotone) metric [H [2T]. The computational challenges involved, 
however, might, at least in certain respects, be even more substantial. 

At this stage of our research, after posting the results above as a preprint, Charles Dunkl 
detailed a computational proposal that he had outlined to us somewhat earlier. The at- 
tractive feature of this proposal would be that it would-holding the exponent k, of |p^^| 

8 



fixed-be able to compute the adjustment factors for general k, rather than having to do 
so for sufficient numbers of individual members of the sequence k = 1, . . . ,N, to be able 
to successfully apply the Mathematica command FindSequenceFunction, as had been our 
strategy beforehand. The proposal of Dunkl (see Appendix) involved parameterizing 4x4 
density matrices in terms of their Cholesky decompositions. The parameters (ten in number 
for the two-rebit case and sixteen for the two-qubit case) would be viewed as points on the 
surface of a unit (due to the trace requirement) 10-sphere or 16-sphere. The squares of the 
points lie in a simplex. One can then employ the corresponding Dirichlet distributions over 
the simplices to determine the corresponding expected values (joint moments). (A further 
facilitating aspect here is that both \p\ and the jacobian for the transformation to Cholesky 
variables are simply monomials.) Using this approach, we were able to extend our single 
{k = 1) two-qubit result ([s]) to the k = 2 case, 

{\pW^?)2-,um/HS = (15) 

fc(fc(fc(fc(fc(fc + 15) + 67)+45) + 220)+4260) + 10944 ^ 
64(2fc + 9)(2fc + ll){Ak + 17)(4A; + 19)(4A; + 21)(4A; + 23) I^-i^^^^/hs- 

Additionally, in the following several arrays, we show (/? = !,..., 12) column-by-column, 

the (3k + 1) coefficients of the numerator polynomials in ascending order-the entries in 

the first row corresponding to the constant terms,. . .-in the two-rebit case. For the cases 
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K= 1,...,6, 

/ -16 4860 -3612816 6610161600 -23680812672000 147885533254368000 \ 
5 2940 -2516616 5496485760 -21644930613600 144374531813568000 



709 


-401334 


1636873812 


-7755993054000 


58524043784903280 


368 


136801 


166748972 


-1199508017652 


11977854861441312 


203 


84291 


6212189 


-4378482660 


1052189083196640 


48 


29493 


13904508 


29246867605 


-30302414250528 


4 


8559 


7805462 


7876634465 


-6899036908859 


— 


1674 


2389416 


2649513956 


3583820785224 


— 


180 


525681 


883461210 


1632448582425 


— 


8 


84496 


219916945 


477741210624 


— 


— 


9112 


40679505 


118164517947 


— 


— 


576 


5660714 


23817008856 






16 


575800 


3786901675 








40000 


469728096 








1680 


44685468 








32 


3143808 










153360 










4608 










64 
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1478171868716632320000 
1513452218833263744000 
-666431742762610272000 
-159042866967850127040 
-20322562269329104680 
-810132874858021428 
132995097356746814 
21392495012420189 
2726092518763299 
1004313607984511 
339532730040063 
85862223778653 
17623087217813 
2989984893665 
414832485207 
46447024562 
4134192972 
285783624 
14766416 
535136 
12096 
128 



22213289955937152264192000 
23581102260458975345971200 
11028809359452120997904640 
2917571120318305111773312 
456603746675135611725072 
36986782703412677565888 

46803055033262612808 
-261889976997036533208 
-9109477993443740087 
4010582339036286008 
1188217816129015268 
313924542711356232 
78095538742712398 
16444648671941448 
2889993079367548 
425775247755632 
52413983854433 
5342837634528 
445432573168 
29856957312 
1568760928 
62140928 
1741568 
30720 
256 
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/ -478312140758259028431974400000 \ 
-522398952485281283922478080000 
-255802623523478150875233408000 
-72876271183514797204357862400 
-12968985411931489708456120320 
-1391118724430074132545941424 
-65460659627903052772432944 
3455511572907702205749432 
652248751783303606044024 
25126439305837999036933 
1467731023076566955385 
1183997600031628820142 
383843332752788856636 
92094572188997438886 
19032858442918448262 
3399681068905043760 
519692040800245908 
67749101821899549 
7505505996250329 
702038674797138 
54914983170480 
3545548122432 
185537763936 
7665748416 
240406272 
5370624 
76032 

V 512 ; 



12 



' 14207139636708483077514789888000000 ^ 
15876402339932769370618323302400000 
8060699595710982594644066482560000 
2429195772959924259673598495040000 
473246195907047274024301397419200 
59861458630180990183509791530560 
4410689910607254768908441542704 
77978523981901390446515581440 
-17699546304068611844956276620 
-1448009072942187597584330860 
23679691892051020702374101 
12657804955771967753837040 
2173418236785245059311455 
535603226600605770029160 
130916289625347750649830 
27174211939515057439440 
4857622320912622046130 
755231957874779338260 
101981982088141404465 
11919937134729947040 
1201281021189026955 
103825857120570960 
7638091703520180 
473427510762240 
24386539757280 
1024463577600 
34174797696 
870092800 
15865600 
184320 

^ 1024 y 
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/ -564438487027912271386331788873728000000 \ 
-642782791049172077649575252852121600000 
-336051652895031776265874186953960960000 
-105851053751774412125536281940889088000 
-22067056295744606471446676944123680000 
-3120323479792255728315846426979831680 
-287292683618793826013052599541333984 
-13768159263526427437128637337477040 
228565888652134052234795741545296 
77614095705972461101015836512688 
3520960278252864643823054097714 
-97356045409654790379642400775 
-3896041249227795478575624541 
3771350531721559074253693217 
1046432219183518743827391855 
230951034077300006360814220 
^ 46540598189766295411264050 / 
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^ 8312454155170970644087050 \ 
1303216963785678085486320 
179484069218701889270145 
21704356449181976554215 
2299019613843658502445 
212567809949303078655 
17076597186282818970 
1184622680244012900 
70388879133182280 
3544333411554144 
149113361559360 
5141635238784 
141440701952 
2982935296 
45253120 
439296 

^ 2048 y 



15 



^ 29240315817392452748492564980040908800000000 \ 
33830315829357716648177866950433062912000000 
18120129848373902187902521023513958348800000 
5915289116881436164243094987941092480000000 
1300158002309897005669893560164992806400000 
199467829529473754135273564780187569164800 
21138278805559013390911138091825622035520 
1413068924132300889731933576834973364800 
35638858522310193185185202018465361456 
-3080272107483898872380758495919446848 
-325601356962501335852286165697351428 
-5934234591428850232759442192641548 
863269405306298197307322635081845 
73444447458140630272343776668948 
7520370028138557798308974457658 
2019677853238030375854720795648 
482836144179289283121850076439 
^ 96138633976278230621072341080 y 
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/ 16899602534589161492520475280 \ 
2646502086250592153436542580 
368329904338448688861367395 
45494778945843226030218660 
4981174263184321953787410 
482464252353850997987160 
41215604000870347678545 
3093092846497823802480 
202871325800172218760 
11552186627193516480 
566265590449861104 
23631870823807488 
827575382720768 
23849804402688 
550461669120 
9776443392 
125356032 
1032192 
4096 



V 



(We are presently attempting extensions to the 13, 14.) The leading (highest- 

order) coefficients in the twelve sets of two-rebit results immediately above are expressible 
in descending order as 



3k+1 



3k 



3 X 2''"'k(k + 2); C3«_i = 2''"^fi:(/t(K(9K + 32) + 24) -45); (22) 



C3«_2 = 2''-^K {k {k [k (9/€^ + 42k + 52) - 119) - 52) - 60) . 



(23) 



From these four formulas, we are able to reconstruct in = 1) all four entries in the first 



column of ( 16 ). Thus, it appears that, in general, Csn-i is a polynomial in k of degree 2(^+1) 



(For i = 3k — 1, we obtain the constant term, of strong interest. With the knowledge of only 
this term, and none of the other coefficients, we can obtain {\p^'^\'^)2-rebit/HS-) Further, we 
have found that 

= (24) 
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K 



1) (135k^ + 855fi:^ + 1895^^ - 1771fi:^ - 3091fi:^ - 7731^^ + 32394^:) 



and 



3k-4 



(25) 



k{k{k{k{k{3k{3k{9k + 59) + 377) - 2887) - 2295) - 10535) + 112240) - 181492) + 436720. 



The numerators of our four sets (k = 1,2,3,4) of two-qubit results are expressible, in 
similar fashion, as 

^ -42 10944 -6929280 9247219200 
-1 
6 
1 



V 



4260 


-3684384 


6039653760 


220 


-456948 


1342859616 


45 


80168 


64072440 


67 


27783 


-13235252 


15 


5373 


1080858 


1 


1458 


1160375 




282 


278478 




27 


50991 




1 


7542 






749 






42 
1 




of all four numerators 



(26) 



character, while the next-to-leading coefficients fit the pattern Cs^ = 3k,{k + 3)/2. 

With our expanded computations-pursuing the Cholesky-decomposition ansatz of DunkI 
we are now able to quite substantially extend the sequence Q to 

((|P||P^^I)'=)2 



0, 



-rebit/HS _ 
{\p\^^)2-rebit/HS 

77 24 209 598 3929 8432 9513091 193880 23471937 1880 22056540991 



(27) 



54' 55' 175' 833' 3724' 9867' 9555975' 211497' 24088922' 1989' 2276223313 J 
{0, 1.4259, 0.43636, 1.1943, 0.71789, 1.0550, 0.85457, 0.99551, 0.91670, 0.97439, 0.94520, 0.96900}. 



Further, we now have for the two-qubit analogue of this sequence 

3 31 _839 35^ 
2' 10'~490' 1260 / 



{Mp''^\)'')2-rebit/HS _ { 3 31 839 3559] 



^p\^^)2-quUtlHS 



(28) 
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{-1.5,3.1,-1.71224,2.8246} 



It is also evident at this point, in striking analogy to the general two-rebit formula (10), 
that in the two-qubit scenario, 

(IpI IP T)2-qubit/HS = -f-{\p\ )2-qubit/HS, (29) 

where, again, both the numerator A^^ and the denominator are 3fi;-degree polynomials in 



k, and (cf. ([TT])) 

B, = 2'-(k+l] (2k + ^] . (30) 



2k 



In the course of this work, Charles Dunkl further communicated to us a result (following 
his joint work with K. Zyzckowski reported in [22], where "the machinery for producing 
densities from moments of Pochhammer type" was developed) giving the univariate proba- 
bility distribution over t G [0, 1] that reproduces the Hilbert-Schmidt moments of t = 2^\p\, 
where p is a generic two-rebit density matrix. (It would be interesting to try to extend the 
methodology employed to the two-qubit and other higher-order cases. Dunkl commented 
that "The formula is slightly misleading near t = 1, there the density is (1 — t)2 times an 
analytic function, I imagine a polynomial approximation is better for computation there, 
but it's obviously the stuff near zero that's important.") This probability distribution took 
the form (cf. ^[eq. (4.3)]) 
63 



(^l-Vi (-8t - 9v^ + 2) + IStlog (^\Jl~Vt + 1) - ^^logW) 



(31) 



(see Appendix below for further details). 

Of course, one may also consider issues analogous to those discussed above for bipartite 
quantum systems of higher dimensionality. To begin such a course of analysis, we have 
found for the generic real 6x6 ("rebit-retrit") density matrices (occupying a 20-dimensional 
space) the result 



,k\ PT\x _ 4fc5 + 40fc^ + 95A;3- 220P-1149A;-1170 ^ 

I IP \)rem-retr^t/HS - 576(A; + 4)(3A; + 11)(3A; + 13)(6A; + 23)(6A; + 25) )-''^*— 

(32) 



Increasing the parameter k from 1 to 2, we obtained that the rational function adjustment 
factor for {\p\^\p^'^\'^)rebit-retrit/HS IS the ratio of 

16/c^+336/c^+2616fc^+8496A;^+12069A;^+101979A;^+903539A;='+3316809A;2+5620320A;+3715740 

(33) 
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to another ninth-degree polynomial 



331776(A; + 5)(3A; + ll)(3A; + 13)(3fc + 14)(3A; + 16)(6A; + 23)(6A; + 25)(6A; + 29)(6fc + 31). (34) 

Additionally, for the generic complex 6x6 (qubit-qutrit) density matrices (occupying a 
35-dimensional space), we have obtained the result 

_ + 15A;^ + 37A;3 - 423A;2 - 2558A; - 3840 

(Ipl IP \),um-,utrit/HS - 72(2A; + 13)(3A; + 19)(3A; + 20)(6fc + 37)(6A; + 41) )1-^^^-'i-tr^t/HS■ 

(35) 



Appendix A: Derivation by C. Dunkl of probability distribution (31) over t € [0,1] 
having the moments t = 2^\p\ 



Notes on moments, etc. C. Dunkl 4/11/11 
Cholesky decomposition: 

Let C be a real upper-triangular N x N matrix, entries Cij, Cij = for i > j and cu > 
for all i. Let P = C^C, entries Pij = X^^i CfciCfcj = YlT=i^'''^ CkiCkj- Consider the Jacobian 
matrix ^ where the dependent variables are Pij^i < j. Claim: 

N 

2-11. 



det 



dp 



dc 



c: 



N+l-i 



i=l 



Lemma: suppose i/i = fi {xi,X2, ■ ■ ■ ,Xi), 1 < i < N then the matrix j is lower-triangular 
(Oforj>.)and det (g;) = nf=i ||- 

This applies to Cholesky: order the variables: Cn, Ci2, . . . , Ciat, C22, • • • , C2n, C33, 
. . . cjy^i^N-i, cn-i,n, cnn- For i < j, Pij = Yl'k'=i CkiCkj + CiiCij] the lemma shows 



det 



dp 
dc 



N N r. N 

nn|^=n(2-.r«) 

1=1 j=i ■' 1=1 



Consider random variables, values in < t < 1. Moments for the Beta distribution: let 
a, /3 > then 



5(a,/3) Jo 



ee-'{i-tf'^ dt 



[a] 



(« + /3). 



,n = 0,l,2. 
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Consider D — detP where P is a random positive definite 4x4 matrix, trace 1. We have 
(expectation) for n = 0, 1, 2, . . . 



(i).(2)„(i)Jl)„ ^ 2-^"(4),„2-^"(2),, 
(10)4n 2-(5),„(f),„ 

1 (4)2„(2)2n 



Let X — 2^D; X is (equidistributed as) the product of two independent random variables 
Xi, X2 with 

(4k _ 4 _ 2 



^ (5)2, 4 + 2n 2 + n' 

E (X^) = 



(2)2„ 



(-) ■ 

Clearly has the density fi (t) ^ 2t,0 < t < 1. The density of X2 is 

1 / M 5/2 



1 



B (2, i) Jo 

(2)2. 



s^"s{l-sf^^ ds 



(¥) 



2n 



The density / {t) of X1X2 is given by 



/ it) = /" /i f -1 h (s) 



2B (2,1) 



I) Jt s 



B (2 
63t 



^ fu-^{l-ufdu 
C u-^l-ufdu. 



2 ./V* 
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The integral is evaluated as follows: set u — 1 — s , du — —2sds, 



ry/l-Vi „6 

f{t)^63t ^ds 

Jo U 



s 



Q3t j -s (15 - 25,s- + 8,s^) ^ 15 (1 + sf ^ ' 



s=0 



8 [ 2 l-s2 

= y Hl-\/t) (^2-9\/t-8t) +15nnM + yi-\/tj - — Uni I . 

This can be easily plotted. Also f (t) ^ O (^{1 - tf^^^ near t = 1. 
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